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A propagation of an ultra-intense short laser pulse in a fiber is investigated with two dimensional 
Particle-in-Cell simulations. The fiber is a narrow hollow channel with walls consisting of overdense 
plasma. In the nonlinear interaction of the laser pulse with fiber walls high order harmonics are 
generated. Sufficiently high harmonics, for which the fiber walls are transparent, propagate outwards 
at certain angle. This is a scheme of a generator of ultra-short pulses of coherent light with a very 
short wavelength. 
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Generation of high order harmonics of the electromagnetic radiation during an interaction of an ultra-intense laser 
pulse with underdense and overdense plasmas is a manifestation of one of the most basic nonlinear processes in physics. 
High order optical harmonics have been observed in the laser plasma interaction for the laser intensity ranging from 
a moderate level to relativistic intensities. 

High order harmonics attract a great attention due to a wide range of their applications for the diagnostics, the UV 
and coherent X-ray sources, lithography etc. (see Refs. 0,0, 013)- Recently, the high order harmonics polarization 
properties have been used in Ref. [j| to detect the strongest magnetic field generated in the laser plasmas. 

Different physical mechanisms of high order harmonics generation have much in common because of a property 
of nonlinear systems to react in an anharmonical manner to a periodic driving force. On the other hand, specific 
realization of this property depends on the circumstances of the laser-matter interaction and mainly on the laser 
intensity. At moderate intensities, generation of high order harmonics is caused by an atom anharmonic response to 
an oscillating electric field. At relativistic intensities, the electron quiver energy becomes higher than the rest mass 
energy and the Relativistic Nonlinear Optics comes into play 0]. In this regime high order harmonics generation is 
due to a nonlinear dependence of the particle mass on its momentum and due to electron density modulations in 
the electromagnetic wave. In the present paper we discuss high harmonic generation in collisionless plasmas by a 
relativistically strong electromagnetic wave. 

In the underdense plasmas an ultra-intense laser pulse causes high harmonics generation through a parametric 
excitation of electromagnetic and electrostatic waves with different frequencies. In an interaction with overdense 
plasmas the laser radiation (partially) reflects back at the plasma-vacuum interface, in the case of sharp plasma 
boundary, or at a critical density surface, in the case of gradual density profile. The reflecting layer in the plasma is 
dragged by the electromagnetic wave back and forth in the direction of the incident light, as well as in the transverse 
direction, forming an oscillating mirror (see Refs. 0, 0, 0, 0, El El ^2 ) • The light reflected by the oscillating mirror 
contains odd and even harmonics, whose polarization and amplitudes depend on the incidence angle, intensity and a 
polarization of the laser pulse. 

High efficiency of the laser energy transformation into the energy of high harmonics can be achieved in a scheme, 
in which the laser pulse propagates inside a hollow channel El- In this case, if the channel is sufficiently long, the 
laser pulse undergoes a multiple reflection at the plasma-vacuum interface, and a portion of the laser pulse energy 
transforms into the energy of high harmonics at each reflection event. If the laser pulse propagates inside a hollow 
channel in a tubiform fiber, the radiation with sufficiently high frequency propagates through the fiber walls. The 
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fiber walls are transparent for harmonics, whose indices are greater than some critical value. Only these harmonics 
are emitted from the fiber at certain angle with respect to the fiber axis. This approach allows one to design a high 
frequency radiation source with controlled properties, which depend on the fiber diameter and the fiber wall thickness. 
We call this source "the relativistic whistle" . 

In this paper we investigate a propagation of an ultra-intense short laser pulse in a fiber, with an aim to study the 
properties of the high harmonic radiation. Here we use the two dimensional approximation, leaving the investigation 
of the 3D problem for forthcoming publications. 

II. MATCHING CONDITIONS 

In the two-dimensional case the fiber corresponds to two thin foils of a finite length. The foils are parallel to each 
other, and to the x-axis in the x, y-plane. Let the wall thickness of the foil is I, and a distance between the foils equal 
to L. We consider the foils to be made of the collisionless overdense plasma with the electron density n e . Such a 
configuration can be formed by an ultra-intense femtosecond laser pulse interacting with the solid density fiber, when 
the fiber wall material is ionized during the first half-period of the laser pulse, and the length of the laser pulse is 
much shorter than the timescale of a hydrodynamic expansion of the plasma. In the framework of this approximation 
we can assume the ions to be at rest. 

From the dispersion equation, to 2 = k 2 c 2 , where u> is the wave frequency and k 2 = k 2 + k 2 , we find the group 
velocity of the wave inside the channel, v gr — dto/dk x = (c/uj)(uj 2 — n 2 /L 2 ) 1 / 2 = ccos6>, where tan 6* = k y /k x = 
(4L 2 /A 2 — l) -1 / 2 and A = 2-kc/lo. The matching conditions for the wave inside the channel, and the wave radiated 
outwards, show that all the outgoing waves propagate at the angle 9 with respect to the channel axis, independently 
of the wave frequency. 

The fact that all the outgoing waves have the same angle of propagation can be demonstrated in a different way. 
After the Lorentz transform into the reference frame moving at a speed of v gr alongthe a;— axis, the a;— component 
of the fundamental mode wave vector vanishes: k' x = (see for details Refs. 0, 13, El)- The problem becomes 
one-dimensional. Therefore high frequency waves excited by the fundamental mode at the plasma-vacuum interface 
also have zero x~ component of the wave vector in this reference frame. They propagate outwards in the direction 
perpendicular to the fiber walls. Performing the inverse transform back to the laboratory reference frame, we obtain 
that all the outgoing waves propagate at the same angle 9 with respect to the channel axis, independently of the wave 
frequency. Their energy is localized inside a slab co-moving with the laser pulse along the x— axis. The slab length is 
equal to the length of the laser pulse. Its transverse size depends on the effective length of the interaction. According 
to the selection rules of the harmonic generation at each wall (see Refs. the s— polarized fundamental mode 

in the fiber generates s— polarized odd harmonics and p— polarized even harmonics of radiation. The p— polarized 
fundamental mode generates only p— polarized harmonics of radiation. 

The amplitude of the n-th harmonic can be found in the approximation of an infinitely thin walls of the fibre 
(see Ref. g). It is E^ = K^E^, where E^ and £W are amplitudes of the fundamental mode and of the 
n-th harmonic, respectively. The transformation coefficient K^ 0,n ' is given by K^ '^ w [(ao/e) sin6>]", in the limit 
nio/ujp e <C e. Here ao = eE^ /m e uJoc is the dimensionless amplitude of the fundamental mode, s = 27rnoe 2 1 /m e uJoc 
is the dimensionless parameter (relativistic transparency measure), too is the frequency of the fundamental mode. 

III. PARAMETERS OF THE COMPUTER SIMULATIONS 

We perform computer simulations using the two-dimensional version of Particle-in-Cell relativistic electromagnetic 
code REMP, based on 'density decomposition' scheme [l5|. Computation box has a size 45A x 20A. Since the problem 
of high harmonic generation requires substantially high spatial resolution, the computation mesh has 50 cells per A. 
The number of quasiparticles per cell is equal to 16 in the plasma region. The channel walls have a length l c h = 25A 
and width 2A, a distance between the walls is L — 0.7A. The ions are assumed to be immobile and the electron density 
is equal to 8n cr , where the critical density is n cr — Wgm e /47re 2 . Since the channel diameter is less than the laser light 
wavelength, only the s-polarized electromagnetic wave (whose electric field is directed along the z-axis) can propagate 
inside the channel. For the parameters chosen, only harmonic, whose index is higher than 4, can propagate through 
the walls. This conclusion follows from the boundary conditions for the electromagnetic wave, which is obliquely 
incident on the plasma-vacuum interface, at the angle 9. The s-polarized laser pulse is initialized in a vacuum region 
at the left hand side from the channel entry. The pulse aperture is f/1, its length is equal to 7.5A, and its dimensionless 
amplitude at the 1A focus spot is ao = 3, which corresponds to the intensity / w 10 19 W/cm 2 . These parameters arc 
close to the parameters of the laser in the CUOS at the University of Michigan Q. We notice here that the pulse 
focusing into the 1A spot is discussed in Ref. Q. In our simulation model the entry of the channel is smoothed to 
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FIG. 1: The electromagnetic energy density and the plasma density distribution at t — 15 x 2tt/uio when the pulse enters into 
the fiber. 
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FIG. 2: The electromagnetic energy density distribution in the x,y plane at t — 20 x 2n /u in frame (a), t = 35 x 2n/u in 
frame (b), and at t — 50 x 2n/u in frame (c). 

decrease the laser pulse reflection. Fig. shows the focusing laser pulse and the fiber. We see that the smoothing of 
the channel entry provides almost reflectionless matching of the laser pulse and the channel. 

IV. RESULTS OF PIC SIMULATIONS 

When the laser pulse propagates inside the channel a portion of its energy is reflected back due to the process 
similar to the stimulated Raman scattering in the plasma of the channel walls; a portion of its energy is transformed 
into the energy of high harmonic waves radiated outwards; the remaining part of the laser pulse propagates through 
the channel toward the end of the fiber, then it goes out as a diverging electromagnetic wave. These successive phases 
are seen in Fig. [21 where we present the electromagnetic energy density distribution in the x, y— plane at different 
times. 

In Fig. |3| we show the distribution of the z-component of the electric field in the x, y-plane. We see the laser pulse 
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FIG. 3: The distribution of the 2-component of the electric field in the a;,t/-plane at t = 45 x 2n/uio. 




FIG. 4: Spectrum of the s-polarized harmonics. The dashed line denotes the dependency k y ~ 0.8^. Each peak corresponds 
to a harmonic with index n = (fc^ + fc^) 1//2 . 



inside the channel and the high harmonics aside the channel. In this figure all the harmonics have odd indices and 
they are s-polarized, as noticed above. The spatial region filled by high harmonics has a length equal to the laser 
pulse length and the transverse size (along the y-axis) equal to 5tc sin 6, where sinfl « 1/2 and St « 20 x 2n /uj n . 

We perform the spatial Fourier transform of the z-component of the electric field inside a sub-domain (18A < x < 
28A; 3A < y < 10A) at t = 30 x 2tt/uj q . The resulting spectrum of the s-polarized light inside this sub-domain is 
presented in Fig. 01 We see the odd harmonics with the harmonic index n = (A; 2 + ky) 1 / 2 > 5. The spectrum 
peaks are arranged along the straight line given by k y rs 0.8k x , i. e. the angle of propagation is smaller than 
arccos(l — 7r 2 c 2 /lo^L 2 ) 1 / 2 . It means that the effective channel width is larger than 0.7A. The increase of the channel 
width is apparently due to the action of the ponderomotive force on the electrons at the plasma- vacuum interface. 

From the expression m e c 2 V(l + a 2 ) 1 / 2 ~ 47moe 2 (5y, which corresponds to the balance between the ponderomotive 
force and the force due to electric charge separation, we find that the channel width becomes larger by 2Sy « 0.1A. 
Assuming now the channel width to be equal to 0.8A we obtain the propagation angle 9 e ff s» arctan(0.8) w 38.7°, in 
agreement with the dependence seen in Fig. 

In Fig. [SI we present the distribution of the ^-component of the magnetic field at the same point in time as Fig. [3] 
We see the high harmonics with even indices. They are p-polarized. As noticed above, the fundamental mode of the 
laser pulse inside the channel does not contain the p-polarized component and we see only the high harmonics inside 
and aside the channel, as well as the low frequency surface mode in the vicinity of the outer plasma- vacuum interfaces 
on both sides of the fiber. The even-index high harmonics outside the channel are localized in the same region as the 
odd-index harmonics. 

The spectrum of the even-index high harmonics calculated for the z-component of the magnetic field in the same 
sub-domain and at the same point in time as in the case of the odd-index harmonics is presented in Fig. [S] We see 
the even harmonics with the harmonic index n > 6. As in the previous figure, the spectrum peaks are arranged along 
the same line k y w 0.8fca;. 
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V. CONCLUSION 



The high harmonic generation in the interaction of the ultra-intense short laser pulse with a narrow fiber is demon- 
strated with PIC simulation. Due to nonlinear interaction of the laser pulse with the fiber walls, the quivering back 
and forth electrons form the oscillating mirrors. The high order harmonics are generated as a result of a multiple 
reflection of the laser light from the oscillating mirrors. The fiber walls are opaque for the electromagnetic wave with 
a frequency below some critical value, which is defined by the electron density of the fiber walls and the channel 
width. As a result, only sufficiently high harmonics, for which the fiber walls are transparent, propagate outwards 
at certain angle, independently of the harmonic index. In addition, we observe the excitation of the p-polarized low 
frequency mode, which is localized at the outer surface of the fiber. The longitudinal scale of the low frequency mode 
is of the order of the laser pulse length. It provides an example of the nonlinear rectification of the light. The high 
harmonics energy is localized in the region with a length (along the x-axis) equal to the laser pulse length, and the 
width (along the y-axis) equal to l c htan9 e ff. This region is co-moving with the fundamental mode inside the fibre. 
The amplitude of the n-th harmonic can be estimated as (oo sin9 e ff) n (uJo/LUp e ) 1+2n , where Lu pe = (47rn e 2 /m e ) 1 / 2 . 
For a = 3, (w /^ P e) 2 = 1/8 and 6 eff w 38°, we obtain £(")/£(<>) « 0.35(0.23)™. This proves an effective mechanism 
for generation of ultra-short pulses of coherent light with a very short wavelength. The corresponding device can be 
called "the relativistic whistle" . 



[1] N. Bloembergen, Nonlinear Optics (Mass.: Reading, 1965); Y. R. Shen, The Principles of Nonlinear Optics (Wiley, New 

York, 1984); S. A. Akhmanov, R. V. Khokhlov, Problems of Nonlinear Optics (VINITI, Moscow, 1965); R. W. Boyd, 

Nonlinear Optics (Academic, San Diego, 1992). 
[2] J. Zhou, J. Peatross, M. M. Murnane, H. C. Kapteyn, I. P. Christov, Phys. Rev. Lett. 76, 752 (1996); Y. Kobayashi, T. 

Sekikawa, Y. Nabekawa, S. Watanabe, Opt. Lett, bf 23,64 (1998); P. Villoresi, et al., Phys. Rev. Lett. 85, 2494 (2000); C. 

Altucci, et al, Phys. Rev. A 61, 021801(R) (1999). 
[3] S. V. Bulanov, F. Califano, G. I. Dudnikova, et al., in: Reviews of Plasma Physics (Kluwer Academic. Plenum Publishers, 



6 



New York, 2001), Ed. by V. D. Shafranov, Vol. 22, p. 227. 
[4] G. Mourou, Z. Chang, A. Maksimchuk, et al., Plasma Phys. Rep. 28, 12 (2002). 

[5] M. Tatarakis, I. Watts, E. N. Beg, et al., Nature 415, 280 (2002); M. Tatarakis, A. Gopal, I. Watts, et al., Phys. Plasmas 
9, 2244 (2002). 

[6] S. V. Bulanov, N. M. Naumova, F. Pegoraro, Phys. Plasmas 1, 745 (1994). 
[7] P. Gibbon, Phys. Rev. Lett. 76, 50 (1996). 

[8] R. Lichters, J. Meyer-ter-Vehn, A. M. Pukhov, Phys. Plasmas 3, 3425 (1996). 
[9] A. V. Vshivkov, N. M. Naumova, F. Pegoraro, S. V. Bulanov, Phys. Plasmas 5, 2752 (1998). 
[10] A. Macchi, F. Cornolti, F. Pegoraro, T. V. Liseikina, H. Ruhl, V. A. Vshivkov, Phys. Rev. Lett. 87, 205004 (2001); A. 

Macchi, F. Cornolti, F. Pegoraro, Phys. Plasmas 9, 1704 (2002). 
[11] M. Zepf, G. D. Tsakiris, G. Pretzler, I. Watts, D. M. Chambers, P. A. Norreys, U. Andiel, A. E. Dangor, K. Eidmann, C. 

Gahn, A. Machacek, J. S. Wark, K. Witte, Phys. Rev. E 58, 5253 (1998). 
[12] A. Tarasevitch, A. Orisch, D. von der Linde, P. Balcou, G. Rey, J.-P. Chambaret, U. Teubner, D. Kloplel, W. Theobald, 

Phys. Rev. E 62, 023816-1(2000); I. Watts, et al., Phys. Rev. Lett. 88, 155001 (2002). 
[13] S. V. Bulanov, F. F. Kamenets, A. M. Pukhov, F. Pegoraro, Physics Letters A 195, 84 (1994). 
[14] A. Bourdier, Phys. Fluids 26, 1804 (1983). 
[15] T. Zh. Esirkepov, Comput. Phys. Comm. 135, 144 (2001). 



